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Abstract: 



^ ! We rewrite the loop equations of the hermitian matrix model, in a way which 

. involves no derivative with respect to the potential, we compute all the correlation 

^ ■ functions, to all orders in the topological expansion, as residues on an hyperel- 

^ ■ liptical curve. Those residues, can be represented as Feynmann graphs of a cubic field 

I theory on the curve. 

iu: 



1 Introduction 

We consider the formal hermitian matrix integral: 

Hn 

where M is a x hermitian matrix, dM is the product of Lebesgue measures of all 
real components of M. V{x) is a polynomial of degree d + 1 > 2 called the potential. 
Z is called the partition function. 

Our goal is to compute the large limit, as well as the full l/N"^ expansion, of the 
following formal expectation values (<> is the average computed with the probability 
measure ^q-^''^^^'^ dM): 

Wk(xi,...,Xk) ■.= N''-^( ti ^tr tr ^\ (1.2) 

^ ' ' ^ ^ xi-M X2-M Xk-M/^ ^ ' 
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Wk{xi, . . . ,Xk) := lim Wk{xi, . . . ,Xk) (1.3) 

A''^oo 

where the subscript c means connected part or cumulant. 

When Z is considered as a formal generating function, it is well known [TT that 
the correlation function Wk{xi, . . . ,Xk) has a expansion, also called topological 

expansion, noted: 

oo 

Wk{xi, ...,Xk):=J2 N-^''Wt\xi, ...,Xk) (1.4) 

h=0 

Let us emphasize that in general Z is not a convergent integral, the partition func- 
tion as well as the W^s are to be understood as formal series in the coefficients of the 
potential, see [TT] for details. 

In that formal sense, the expectation value of a product of k traces is the com- 
binatoric generating function for enumerating discrete surfaces with k holes, and the 
variables fugacities for the lengths of the k boundaries [23 El 

It is well known ^T] that the power of associated to each discrete surface is 
its Euler characteristic x- Fo^' ^ genus h connected surface with k holes, we have 
X = 2 — A; — 2h. This is why the quantitv 11.21 has a large limit. The large N limit 
^^Wk{xi, . . . ,Xk) = w'^\xi^ . . . ,Xk) is therefore the generating function of genus 
zero discrete surfaces with k boundaries, and each wl!'\xi, . . . ,Xk) is the generating 
function of genus h discrete surfaces with k boundaries. 

The problem of computing the WkS has been addressed many times, for various 
applications to physics and mathematics. Indeed, the correlation functions of eigen- 
values (and thus of traces of powers) of a random matrix have a universal behavior 
which (this is what universality means) is observed in many physical phenomena, rang- 
ing from solid state physics (quantum chaos, mesoscopic conductors, see [HHl 1211 ) to 
high energy physics (nuclear physics [HOI, Quantum chromo dynamics [321 , string theory 
|12j). and in mathematics (distribution of Riemann zeta's zeroes [HI). 

In the 90 's, random matrices were extensively studied in the context of quantum 
gravity (see [H]), which is nothing but statistical physics on random discretized sur- 
faces, i.e. the combinatorial problem of enumerating discrete surfaces of given topology, 
as described above. 

Quantum gravity is also deeply related to conformal field theory (CFT), when one 
takes a "double scaling limit" where very large discrete surfaces are dominant, in other 
words, CFT is the limit of continuous surfaces. Depending on the limit chosen, and 
on the coefficients of the potential V{x), one may reach different double scaling limits, 
which are in relationships with the {p, q = 2) minimal models in CFT. All the critical 
exponents of such surfaces are given by KPZ's formula j27j . 
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It is thus expected, that in appropriate double scahng hmits, expectation values of 
the form II .21 can be computed from a quantum field theory, namely Liouville's theory. 

Here, without taking any double scaling limit, we will find a quantum-field-theory- 
like Feynmann expansion for the Wk{xi, . . . , Xk)- 

The WkS have been computed in the literature by various methods. A formula of 
Dyson J3] gives the W^s for finite in terms of orthogonal polynomials, but is not 
very convenient for large N limit calculations, and is not convenient for the formal 
model. The method of loop equations [HI EH El IHl UHl HE] gives recursion relations 
between the WkS, which simplify in the large limit. The loop equations have been 
known for a while, and give a very effective algorithm for computing explicitly the 
WkS (see 12 HI El)- The method developed by |Tj for computing the W^s, consists 
in computing Wi and then obtain the W^s by taking iterated derivatives with respect 
to the potentials (loop insertion method). This method has two drawbacks: first in 
order to find Wk+i, one should know Wk for all potentials (in particular one must 
take infinite degree potentials); second, before computing Wk, one has to compute 
Wi, W2, . . . , Wk-i, i.e. it has not been found how to integrate the recursion formulae 
of p. 

Here we consider new loop equations, which allow to find recursion relations between 
the WkS, without taking any derivatives with respect to the potential (we may work 
with fixed potential). Moreover, the recursion relations for the wj^^^^s obtained in 
this paper can be integrated: the k^^-loop function to order A^"^'^ is a /c-legs, /i-loops 
Feynmann graph of a 0^ theory living on an hyperelliptical curve. 

Outline: 

- in section 121 we introduce the notations. 

- in section jSl we introduce some basic tools of algebraic geometry. 

- in section |3] we write and solve the loop equations to large A^ leading order, i.e. 
we compute the W^s. 

- in section we write and solve the loop equations recursively to each order in 
l/N"^, i.e. we compute the W^'^'^^s. 

- in section ini we do explicitly the computation in the one-cut case. 

- in section Owe conclude by presenting perspectives of applications to other matrix 
models (2 matrix model). 

2 Definitions and notations 

From now on, we assume that V'{x) is monic of degree d > 1. 
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2.1 Loop functions 

For /c > 1, we define (the subscript c means connected part or cumulant): 

Wk(xi,...,Xk) ■.= N^-UtT ^tr tr ^\ (2.1) 

rk-2/.VIix,)-VliM,)^^ 1 1 



UJxi; xs, . . . , Xk) := N"-' ( tr '-^^ tr . . . tr ) (2.2) 

^ ' ' ' ^ \ xi-M X2-M Xk-M/^ ^ ' 

and their large N limits: 

Wk{xi,...,Xk) ■.= Wl°\xi,...,Xk) := Mm Wk{xi, . . . ,Xk) (2.3) 

Uk{xi;x2, ...,Xk) := U^^\xi;x2, ...,Xk):= lim Uk{xi;X2, ...,Xk) (2.4) 
as well as their formal expansions {k > 1, h > 0): 

oo 

Wk{xi, ...,Xk):=J2 ■■■,Xk) (2.5) 



h=0 



Uk{xu X2, . . . , Xk) := N-^^ Uj!'\x,; X2, . . . , Xk) (2.6) 



h=0 



Notice f/f ^ is a monic polynomials of degree d — 1, and as soon as k + h>2, f/f ) 
is a polynomial of degree at most c? — 2 in xi. We have: 



[/f)(xi;x2,...,Xfc) = Pol Vl{xi)wi^\x,,...,Xk) (2.7) 



The functions VTfc are called loop-functions, because they are generating functions 
for discrete surfaces with k boundaries, i.e. k loops. 

2.2 Filling fractions 

If the integral ll. II were to be considered as a convergent integral, the expansion 
would exist only in the so-called one-cut case (see [71[TU]). Here 11.1 1 is considered as 
a formal power series, by its expansion in the vicinity of a minimum of the potential 
tTV{M). The potential V{x) has in general d = deg V extrema, and thus, the potential 
Tr V{M) can have extrema indexed by the number of eigenvalues of M lying in the 
vicinity of each extrema of V{x). The formal perturbative expansion around such local 
extrema cannot change the number of eigenvalues near each extrema. The fractions 
of eigenvalues near each extrema of V are called filling fractions, and are thus moduli 
characterizing the vacuum near which the perturbative formal expansion is computed. 
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The filling fractions play an important role in recent applications of random matrix 
models to string theory |12j . 

The filling fractions are denoted: 

s 

ei,e2,...,e, , ^ej = l (2.8) 

i=i 

It is well known (and we recover it below) that the function Wi{x) is solution 

of an algebraic equation, it has s cuts [a2j_i, a2j], j = 1, . . . , s, which correspond to the 
location of eigenvalues in the large limit. The condition that the filling fractions are 
given can be written: 

Vj = l,...,s, Wi{x)dx = e, (2.9) 

where the contour surrounds the segment [a2j-i,(i2j] in the trigonometric direction. 

Let us for a moment, use the method of for finding the filling fraction condi- 
tions for other loop functions. 

If V{x) = 'Ylik^kX^ ) from |1| we introduce the loop insertion operator^: 

d ^ I d 



dV{x) ' ^ x^+i dt 



we then have P): 

d 

Wk+iixi, ...,Xk, Xk+i) = TTTTT r Wkixi, ...,Xk) (2.11) 

dV{xk+i) 

and thus, since the filling fractions are given parameters independent of V and A^, we 
must have for all > 1, /i > 0, + A; > 1: 

Vj = l,...,s, i" Wl''\xi,x2,...,xk)dxi = (2.12) 

J [a2j-l,a.2j] 

From the same argument, since we assume that there is no eigenvalue elsewhere in the 
complex plane, we can write, for any m in the complex plane, away from the cuts: 



iyi(a;i)da;i = (2.13) 
(where the contour integral is a small circle around m) and thus: 

/ VrfHxi,X2,...,Xfc)dxi = (2.14) 

J m 

^The loop insertion operator j-fc+i is a formal notation which makes sense order by order 

in the 1/x expansion, and eg 12. 121 is perfectly rigorously proven. 



5 



3 The one-loop function and algebraic geometry 

3.1 The one— loop function 

It is well known ^T] |2] (and it is re-derived below) that the one loop function is 
algebraic: 

Wi{x) = ^ {y'{x) - M(x) (3.1) 

U^[x) = i {V'\x) - M\x)(t{x)) (3.2) 
where M and a are monic polynomials (remember V'{x) is monic), determined by: 

Wi{x) ~ - (3.3) 

x^oo X 

and bv l2.9[ which can be rewritten as follows: let ai, . . . , a2s be the zeroes of cr: 

2s 

a{x) = Y[{x - ai) (3.4) 

1=1 

we must have: 

Vj G [1, s - 1], / M{x)^/^ix)dx = 2mej (3.5) 

For a given s G [1, ci], the eg nations 13 . 31 and 13 . 51 give a finite number of solutions for M 
and a. Let us assume that we have chosen one of them. 

3.2 More notations 

For convenience we introduce mi, ... , nid-s the zeroes of M: 

d—s 

M{x) = Y[{x - mi) (3.6) 

i=\ 

We also define for /c > 1, /i > 0, and /i + /c > 1: 

k 



Ff(xi,...,x,):=(2'^iy,(xi,...,a;,)+ ^^'^'^^^'^ 



(xi - X2)'' 



) n (3-7) 



and: 

.(0) 



Ffc(xi, . . . ,Xfc) := \xi, . . . ,Xfc) (3.8) 

It is well known that the -F^'s and Fjf^^^s are rational functions of all their arguments 
(see P 121 El). 
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Another useful notation is in terms of multi-linear differential forms: 

Gk{xi, . . . , Xfc) := Wk{xi, . . . ,Xk) dxi. . . dxk (3.9) 

and for higher orders: 

Gf\xi, ...,Xk) := W^^\xi, . . . dxi ...dxk (3.10) 

It is well known that they are all multi-linear differentials defined on an hyperelliptical 
surface. All of them, except Gi and G2-, have poles only at the branch-points (i.e. the 
zeroes of a), and have vanishing contour integrals around the cuts. 
All this is re-derived below. 



3.3 Hyperelliptical surfaces 

We need to introduce some basic notions of algebraic geometry ■ 

Equation 13.11 defines an hyperelliptical surface of genus s — 1. Let y = V'{x) — 
2Wi{x), we have: 

= M^{x)a{x) (3.11) 

That equation defines a Riemann surface with two sheets (corresponding to the two 
determinations of the square root). In other word, for each x, there are two values of 
y{x). 

Let us define the physical sheet as the sheet where: 

x-'^/^ ~ +1 (3.12) 

and the second sheet as the one where: 

x-'^/oix) ~ -1 (3.13) 

If X is a point in the physical sheet, let us note x the point corresponding to the same 
X in the second sheet. By definition, we have: 

= — a/(t(x) , M(x) = M{x) , y(x) = —y{x) , dx = dx 

(3.14) 

The branch points are the points where the two sheets meet, they are such that: 

Vz = l,...,2s, lii = tti (3.15) 
Near a branch point Oj, the surface is better parameterized by the local coordinate: 

Ti{x) := y/x — Qi = —Ti(x) (3.16) 



7 



I.e. 



X = tti + T, , dx = 2TidTi 



(3.17) 



In particular, the differential dx has a (simple) zero at x = Oj. 

Holomorphic differentials: Let L{x) be any polynomial of degree < s — 2. Since 
\/oTx) has a simple zero at x = a^, the differential Lix)—p= has no singularity on the 



whole surface (neither near the branch points, nor at oo), it is thus called a holomorphic 
differential. One has the following classical theorem: there exist a unique set of s — 1 
polynomials of degree s — 2, which we note Lj{x), such that: 



Wl,je[l,s-lf 



[a2l-l,a.2l] 



(3.18) 



The differentials 



dx are called the normalized holomorphic differentials. Notice 



that the Lj's form a basis of degree < s — 2 polynomials. For any polynomial P{x) 
such that deg P < s — 2, we have: 



P{x') 



dx' I Lj{x) 



(3.19) 



Notice that on the s*^ cut, we have: 
VjG[l,s-l], 

We define: 



[a2s-ua2s] \/oix) 



Lj(x) , 

^ dx = -2i'K 



Ls{x) := 

so that 13.1^ holds also with the sum on j running from 1 to s. 



(3.20) 



(3.21) 



Normalized differential of the third kind: For any x' on the curve, there exists a 
unique meromorphic differential, noted dS{x,x'), which has only two simple poles in 
X, located ai x = x' and x = x' and such that: 



dS{x^ x') 



dx 



x^x' X — X' 

dS{x,x') ~_ 

x—*x' X X 



+ finite 



finite 



(3.22) 



Vj 



s-1. 



dS{x, x') 







[«2j-l,a2j] 



Notice that dS{x, x') is a meromorphic differential in the variable x, and a multi-valued 
function of the variable x' . 
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It is easy to check that we have the following expression: 

s-l 

X — x' 



dS{x,x') = f ^ - ] dx (3.23) 

where 

CJx') := 7^ / (3.24) 

In this formula, it is assumed that x' lies outside the contours [a2j-i, a2j]. One has to 
be careful when x' approaches some branch point Uj. When x' lies inside the contour 
around [a2j-i,a2j], then one has: 

Q{x') + = ^ / d^^ 25) 

which is analytical in when x' approaches a2j-i or a2j- 
For i = 1, . . . , s, we define: 

Lj(x) 



dS2i-i{x,x') := dS2i{x,x') := dS{x,x') 



which is a one-form in x, with poles at a; = a;' and x = x', and which is analytical in x' 
when x' is close to a2i-i or a2j. 

Bergmann kernel: For any x' on the curve, there exists a unique bilinear differential, 
noted B{x. x'). called the Bergmann kernel, which has only one double pole in x, located 
at X — x' (in particular no pole at x = x'), with no residue, and such that: 

B[x, X ) ~ + finite 

[x-x) ^3 27) 

Vj = 1, i B{x,x')^Q 

J x£[a2j-i,a2j] 

It is easy to check that B{x, x') — B{x', x) and: 

B{x,x') = dx dx' A ( + _ g c,ix')L,ix)V^) 



7(^ + ^^(^>^')) (3.28) 



It can be written: 



" 2{x - x'y ^ 4(a; - x')^^/^^/^ 
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where Q{x,x') is a symmetric polynomial in x and x' , of degree at most s, such that 
Q{x,x) = 2a{x) and dx'Q{x,x')\x'=x = (t'{x): 

Qix, x') = 2(t(x) + (x' - x)a\x) + - i^lll^ Six) + 0{x' - xf (3.30) 

3 2 

where ^(a;) is called projective connection at x. We can write: 

Q{x, x') = 2a{x) + {x' - x)a'{x) + {x - x'fA{x, x') (3.31) 

where A{x, x') is a polynomial in both variables. We have: 
B{x,x') 1 ^ cr(x) 



dxdx' 2{x-x'y 2(x-x')2v/^v^^M 



+ (3.32) 



4(x' - x)^/cr{x)^/a{x') 4a/ct(x) a/oOtO 

4 Loop equations 

Now, we will introduce a method for computing the W^'^'^^s. It is based on the so-called 
loop equations or Schwinger- Dyson equations, i.e. invariance of the integral ll. II under 
local infinitesimal change of variable. 

4.1 Useful notations 

Let K = {2, . . . , k}. For any j < A; — 1 we denote: 

K,:={ICK /i^I = j} (4.1) 
and for any subset / G Kj we define: 

as well as: 



v'^:=nV^K) (4.3) 

(=1 

and 



dxr.= Yldx,^ (4.4) 
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4.2 Loop equations 

The invariance of the matrix integral ll . II under the change of variable M — > M + ri5M 
(see im UHl Unjfor detailed derivations): 



= n tr 4.5 

^ i=2 ^ 



implies, to first order in r]: 

k=l: W,{xif + -^W2{x,,xi) = V'{x,)Wi{xi) - f7i(xi) 
k>2: 2Wi{xi)Wk{xi, ...,Xk) + -^Wk+i{xi,xi,X2, ...,Xk) 

fc-2 

+ ^ ^ Wj+i{Xi,Xi)Wk-jiXi,XK-l) 
j=l l€Kj 
k 



+E 



d Wk-l{x2, ...,Xj,...,Xk)- Wk-l{x2, ■ ■ ■ ■ ■ ■ 



J= 

V'{xi)Wk{xi, . . . , Xfc) - Uk{xi; xs, . . . , Xk) (4.6) 



i.e., to leading order in l/N"^ we have: 

k = l: Wi{xif = V\xi)Wi{xi) - Ui{xi) 
k>2: 2Wi{x,)Wk{xi,...,Xk) 

k-2 



+ ^ ^ Wj+i{Xi,Xi)Wk-jiXi,XK-l) 

j=l IdK, 
k 

^ ^ d Wk~l{x2, . . . ,Xj, . . . ,Xk) - Wk^l{x2, . . . ,Xi, . . . ,Xk) 

~. dxi 

= V'{xi)Wk{xi, ...,Xk)- Uk{xi] X2,..., Xk) (4.7) 

Notice that it implies ITT] for k = 1. 

Now assume k >2, and using ITT] we rewrite: 

M{xi)^/a{xi)Wkixi, . . . ,Xfc) 

fc-2 

= ^ ^ Wj+i{Xi,Xi)Wk-j{Xi,XK-l) 
j=l IdK, 

^ d Wfc_i(x2, . . . , Xj, . . . , Xfc) - Wk-l{x2, . . . , Xi, . . . , Xfc) 

^ dxj Xj — xi 

+f/fc(xi;x2, . . . ,Xfc) (4.8) 
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4.3 Case k=2 

For k = 2,01 reads: 



0X2 V X2 - Xi J 

1 d ( M{x2)\/(y{x2)-M{xx)\/(y{xx ) 

2 (9X2 \^ X2- Xx 

2 dx2 I X2 — Xi 

1 ^ ( r^^M{x2) - M{xx 



\/(y{x2) 



(4.9) 



2 (9x2 V X2 — xi 

2 (9X2 V X2 - Xi / 



which can be written: 



g{xx)W2{xx, 



X2) 



1 d ( A/cr(x2) - A/cr(xi) \ i?2(a;i;x2) 



2 9x2 V X2 — Xi j M(xi) 

(4.10) 



where i?2(2;i;x2) is a polynomial in xi of degree at most d — 2. From lTT^ we know 
that the LHS has no pole at the zeroes of M, thus 

R2{xi; X2) = M{xi)P2{xi; X2) (4.11) 

where ^2(2^1; X2) is a polynomial in xi, of degree s — 2. We have: 

Z (7X2 \ X2 — Xi / 

(4.12) 

In terms of the function F2 introduced in 13.81 we have: 

F2(xi,X2) _ 1 ^ ^/oix2} ^ P2{xi;X2) ^^^^^ 



4:^/(T{xl)^(r{x2) 2^(r{xi) dx2 (xi - X2) ^(y{xi) 

which proves that F2 is a rational function of Xi, and by symmetry, it is also a rational 
function of X2. 

Then using ITT^ as well as 13. 191 and I3.24| we find P2: 

^'";^ = 1AV^ f^L^_VQ(x2)LKx0l (4.14) 

A,/^^^ 2 9X2 \{X,-X2) ^ ^ ^ 

12 



then using IH7^ 

F2(Xi,X2) 



1 d dS{xi,X2) B{xi,X2) 1 



4:^/(r{xl)^ya{x2) 2 8x2 dxi 



dxidx2 2 ixi — X2) 



(4.15) 



where we recognize B the Bergmann kernel introduced in 13.271 Finally, we have the 
two-loop function in the form: 




or, using 



F2(Xi,X2) 



Q{Xi,X2) 



(4.16) 



(4.17) 



[xi - X2y 

The result I^!TH1 or 14. 1^ is well known and can be found in many places in the literature 
[T3] . We have just presented one derivation for completeness. Now, let us move to 
higher loop functions. 



Remark: we can write 

o[x\ o'[x\ 
F2 x, xi = 2 . ^ \^ ^ + A{x, xi 

[X — X\Y X — X\ 

where A{x^ X\) is a polynomial in both variables. It implies: 

8a{x) ^ 16a(x)2 ^ Aa{x) 
which is a rational function of x, with double poles at the branch-points. 

4.4 k=3 

Starting from 14.81 for k = 3, i.e. 



(4.18) 



(4.19) 



M{xl)^/a{xl)W3{xl, X2, X3) 



2W2iXi,X2)W2iXi,X3) 

d 1^2(3:2, 3:3) - 1^2(3^1, 2:3) 



dx2 X2 — Xi 

^ d W2{X2,X3) -W2{Xi,X2) 

dxs X3 — Xi 

+ U3{Xi;X2,X3) 



(4.20) 



and using the results for A; = 2, we get: 

F3{xi,X2,X3) F2{xi,X2)F2{xi,X3) 



^/oix2) ^/oix3) 



a(xi)M(xi) a/o^X2) A/a(x^ 
13 



^ 8 d 1^2(3:2, X3) ^ 8 d 1^2(^2, X3) 



M(xi) 9x2 X2 — xi M{xi) dxs X3 — xi 

, SU3{Xi;X2,X3) 

+ M(a;i) ^ ^ 

i.e. F^lxi, X2, X3) is a rational function of xi, and by symmetry, it is a rational function 
of all its arguments. Expression 14.201 clearly shows that F3 has no pole when Xi = X2 
or Xi = X3. Moreover, from l2.lil we know that it has no pole at the zeroes of M. Thus 
the only possible poles of -F3 are at the branch points and at 00. Notice that only the 
first term in the RHS of 14. 2 II has poles at the branch points. 
Before continuing, let us study the case k > 3. 

4.5 k larger or equal to 3 



Now assume k > 3. We start from 14.81 

k-2 



^/a{xi)Wk{xi,XK) = 



Wj+i{xi,xi)Wk-jixi, xk- 



Mixi) 

j=l l€Kj ^ 1^ 



k 



+E 

i=2 



d Wk-iixx) - Wk-i{xi,XK-{i}) 
dxi {xi — Xi)M(xi) 



Uk{Xi]XK) 

M{xi) 



(4.22) 



and we consider separately the terms corresponding to j = 1 and j = — 2 in the RHS, 
and we write J2ieKi ~ Si=2' S^^- 



^/a{xi)Wk{x,,XK) = 2^2^— 

^ (W2{Xi,Xi) + 2(cc,-x,r ) Wk-l{Xi,XK-l) 

A d Wk-ijxK) Uk{xi;xK) , . 

^^dx,{x,-x,)M{x,) ^ M(xi) ^ ■ ^ 

This clearly proves, by induction on k, that for all A; > 3, Fi:{xi,xk) is a rational 
function of Xi (and by symmetry, of all its arguments), with poles only at the branch 
points and at 00. We have just re-derived it in a way different from jTll2]- 

Now, assume k > 3. Consider the Euclidean division of the polynomial 
Uk{xi]X2, . . . ,Xfe) by M{xi): 

Ukixi, xk) = I ^ ^ Pkixi, xk) M{xi) + Qk{xi\ xk) (4.24) 
\/cr[xK) 
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where deg Pk = s — 2 and deg Qk < d — s. 

Thus, we have found that for any A; > 3, we have: 

Fk{xi,XK) - Pk{xi;xK) 1 sr^ sr^ Fj+i{xi,Xi)Fk-^j{xi,XK-i) , \ /a nr,\ 

1^=^= = o ^ \- -^ki^i, Xk) (,4.25) 

^/a{xK) ^j^^j^j^^ a{xi)M{xi)^ya{xK) 

where Rk{xi] xk) is a rational fraction of Xi with no poles at the branch points neither 
at oo (it has poles at the zeroes of M and at the Xj's). 

4.6 Cauchy formula 

Cauchy formula gives: 

Fk{xi,XK) - Pkjxu Xk) _ dx Fkjx, xk) - Pkjx; xk) ^5) 

\/<y{xK) x^x^x-xi ^/a{xK) 

where the integrand has poles only at the branch points. Therefore we may deform 
the integration contour used to compute the residue, into residues at the branch points 
only: 

Fk{xi, Xk) - Pk{xi; xk) = Res {Fk{x, xk) - Pk{x; xk)) (4.27) 

1=1 ^ 

using 1^211 we get the recursion formula for the F^s, for all k > 3: 

F,{x„XK)~Pk{x,;xK) = ^ E Res X: E rfx (4.28) 

2 ^ a{x)M{x){xi - x) 

Pk{xi,XK) which is a polynomial in xi of degree at most s — 2, is computed with 
formula 13.191 



Starting from I2.12| i.e.: 

Fk{xi,XK) 



[a2l-i,a2i] 



dxi = (4.29) 



we have: 

Pk{xi;xK)dxi 

[1^21-1,0-21] 

yyyl dx^ F,^,{x,x,)Fk {x xk ,) 

tr U hi, ho2.-..02A ^-''^ 2(xi - x)M{x)a{x) 

^'"^ " dxi f Fj+i{x,xi)Fk-j{x,XK-i] 



tr U h.^-i,a2,\ J[a2.-.,a2^ 2(xi - x)M{x)a{x) 
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dx 



(4.30) 

Notice that for i = I, the contour of integration of xi encloses the contour of x. We 
may exchange the position of the two contours, by picking a residue at xi = x, thus: 

Pk{xi;xK)dxi 

[(12; -1,121] 

" Fj+i{x,Xi)Fk-j{x,XK-i) dxi 



k-2 



j=l iGKj 

(4.31) 



^ ^ ^ / Fj+i{x,xi)Fk-j{x,XK-i) ^^ 



[a2i-i,a2i] 2M{x)^/cr{x)a{x) 



Using the function C/(x) introduced in 13.241 we have: 

Pk{xi;xK)dxi 

Ia2i-i,a2i] y^alxij 

(4.32) 

and thus we have computed Pk'- 



^ ^ " Fj+i(x,xj)Ffc_j(x,xx-j) 



Pk{Xi]XK) = EEy] r OAJt \ t \ 



(4.33) 

4.7 The recursion relation 

That gives the recursion relation for the FkS: 



^ s k-2 

^k{xi;xK) = — J2Y.Y. 



2m 

i=l j=l Ii^Kj 



(4.34) 
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where it is important to remember that the term inside the bracket is analytical when 
X approaches a2i-i or a2i. This allows to write the contour integrals as the sum of two 
residues around a2j-i and a2i. 

It is interesting and more intrinsic to rewrite 14.341 in terms of multi-linear differen- 
tials Gk on the hyperelliptical curve. 

First, notice that a contour around in the hyperelliptical curve (Tj(x) = ^Jx — ai 
around 0) is twice the contour in the complex plane [x around aj), i.e. we will have an 
extra factor 2 in the denominator. 

Then, notice that G2 and F2 differ by a term, which is an even function of the local 
parameter Ti[x\ i.e. which does not contribute to residues near the branch-points (this 
can be checked separately for = 3 and A; > 3). 

Thus we get the recursion relation for the G^'s: 



ri I \ sr^S^ST^ o Gj^x[x,xi)Gk^j[x,XK-i) -ici I 
Gk{xx\ xk) = } y y ^ Res n ^ \ J dS,{xi,x) 

(4.35) 

where now the residues are computed on the hyperelliptical surface (i.e. extra factor 2 in 
the denominator), and dSi{x, x') is the abelian differential of the third kind introduced 
inHSni 

That recursion relation allows to compute Wk in a tree-like recursion from residues 
of lower loop-functions. 

4.8 Solution of the recursion relation as cubic- Fey nmann trees 

Equation 14.351 is conveniently represented with diagrams "a la Feynmann". 

Let us represent the fc-loop correlation function Gk{ black disk with 

k legs: 



:= Gk{xi,. . . ,Xk) (4.36) 




and introduce the following Feynmann rules: 
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Arrowed propagator: 


— X2 := dS^{x,x') 


Vertex: 


( X ■=\l{2y{x)dx) 


Non-arrowed propagator 
=2-loop correlator: 


^1 ^ := " =G'2(a;i,X2) 



Then 14.351 can be represented as follows 




Whose solution is clearly that the fc-loop correlator G^. is the sum over all plane 
binary trees with 1 root and k — 1 leaves, with a skeleton made of oriented arrows 
(oriented from root toward leaves), and whose k — 1 leaves are non-arrowed propagators 
finishing at the x/s with 2 < j < k. 

Notice that two trees which differ only by the ordering of branches at a vertex give 
the same contribution to Gk, so that instead of summing over plane trees, one can sum 
over non-plane trees, with a factor 2*^"^. 

Let 7fc be the set of plane rooted binary trees with k — 1 labeled leaves (0:2, . . . , Xk). 
and let be the set of non-plane rooted binary trees with k — 1 labeled leaves 
{x2, . . . , Xk). We have: 

iVfc+2 := CardTfc+2 = A; + 1! = — = 2^' {2k - 1)!! (4.37) 

kl 

where is the Catalan number which enumerates plane trees. And: 

Oh} 

Nk+2 ■■= CardTfc+2 = 2-'=CardV2 = 2-'— = {2k - 1)!! (4.38) 

kl 

For any given tree T G 7^, with root Xi, leaves xj {j = 2, . . . ,k), and with k — 2 
vertices noted x'^ {v = 1, . . . , k — 2) , so that its inner edges are of the form vi V2 and 
its outer edges are of the form w — j, we define the weight of T as: 

W dSi^{x'^,xJ W G2«,Xj) 

inner edges t;— >ui outer edges v^j 
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(4.39) 



Thus we have: 



(4.40) 



G3 is thus given by A^3 = 1 tree, G4, is the sum of A^4 = 3 diagrams, G5 is the sum 
of iVs = 15 diagrams,... 



4.9 Example: 3-point function 

As an example, let us carry out explicitly the computation for the 3-point function. 
Diagrammatically, we have: 




= 2 




Eg J4. 281 can be written for k = 3: 

F3(xi,X2,X3) = Resdx ^ + -P3(xi, X2, X3) (4.41) 



i=l 



M{x)a{x)ci'{x) {x — x\) 



where P3 is a polynomial in x\. Using OHl notice that F^iai.x) = + A{ai, x) is 
finite, and that: 



I.e. 



a'(x) a(x) ^ , 

^ ' 2- — ^-^-F2{x,xi)+A{x,xi) 



X — Xi {x — Xi)^ 



(4.42) 



J^3{Xi,X2,X3) - 2 2^Kesdx ^^(^.y^._^^^2^;(^ 

i=l 
2s 

— > Res dx 



1=1 
2s 



+ Res dx 



i=l 



M{x){x — Xi)^cr'(x) 

^2(3;, Xi)F2{x, X2)F2{X, X3) 

M{x)a{x)a'{x) 

A{X, Xi)F2{x, X2)F2{X, X3) 

M{x)a{x)a'{x) 



(4.43) 



i?3(Xi,X2,X3) 
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The first line has no residue at the branch-points (indeed, using 14.181 notice that 
F2{ai, x) is finite), and the last line is a polynomial in xi (indeed A{x, xi) is), therefore: 



2s 



N ^ F2{x,Xi)F2{x,X2)F2{x,X3) 
^3{Xl,X2,X3) = - > Resdx . , \ ,/ \ \- ^3{Xl,X2,X3) 



1=1 

2s 



M{x)a{x)a'{x) 



EF2{ai, Xi)F2{ai, X2)F2{ai, xs) 
N T ^3[Xl7 ^2, X3) 



(4.44) 



1=1 



M{ai)a'{a 



where i?3(xi, 0:2, ^3) is a polynomial in xi, of degree at most s — 2. Condition 12.121 
implies that i?3(xi, X2, Xa) = 0, thus: 



F3{xi,X2,X3) 



2s 

E 

i=l 



F2{ai, xi)F2{ai, ^2)^2(0^, X3) 
M(a,)a'(a,)2 



(4.45) 



This is a generalization of what was found in |T1 |2l IHj for the one-cut case s = 1 . 
Let us redo this computation in a more intrinsic way. Start from 14.331 for A; = 3: 



2s 



G'3(xi,X2,a;3) 



2y Res ^^^^'^^^^^f'^^^ dS,{x„x) 



(4.46) 



i=l 



2y{x)dx 



Notice that both dSi{xi, x) and y{x) have a simple zero at x = thus, 
dSi{xi,x) dxdSi{xi,x) 



y{x) 



dy{x 

B{xi, x) + B{xi,x) 

dy{x) 
G2ixi,x 



+ 0{y/x - tti) 



+ 0{y/x — tti 



dy{x) 



+ 0{y/x - tti) 



(4.47) 



This implies: 



2s 



G3{xi,X2,X3) = 2>^Res 

^ — ^ ai 



i=l 



G2{X, X2)G2{X, X3)G2{X, Xi) 

dx dy 



(4.48) 



This agrees with [21113] (our dy is half the dy of [5j). One can also write: 

W3(Xi,X2,X3) = ResW2{x,X2)W2{x,X:i)W2{x,Xl)—— 
^-^ ai rill 



1=1 



dy 



(4.49) 
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4.10 Example: 4 point function 



Diagrammatically, we have: 



'3 = 4 1 




+ 4 X 




+4 




Explicit computation of 14 .341 for k = 4 gives: 



2s 



Fi{xi,X2,X3,Xi) 



F2{x3,aj)F2{x4,aj) ^2(^2, ai)F2(a;i, a^) 



M{aj)a'{a,) 



2s 



Y- F2ixs,aj)F2ix2,aj) ^ 
- 2^ — — TTT — ^2{ai,a_ 



.^.^^ M(a,)a'(a,) 

2s 
2s 



M(a,)(T'(a,)2 

F2(x4,aj)F2(xi,aj) 
^'^^^ M(a.)a'(a,)2 

F2(x3,ai)F2(xi,ai) 



i=l 
2s 



F2(X2, aj)F2(x4, 

M(a,)a'(ai)2 
(F2(x3, ai)F2(a;4, 0^)^2(0^, X2)F2{ai, xi))' 



i=l 



F2{x3, ai)F2{xi, ai)F2{ai, X2)F2{ai, xi) 
M(a,)V'(a,)3 
M'{ai) _ A{ai,a,) 5a"{a,) 
M{ai) a'{ai) Qa'{a,) 



(4.50) 



5 Higher genus 

Now, we don't drop the term in 14 .61 and we expand to order h: 

2f2wt"'\x,)wt\xu...,x,) 



thus: 



m=0 



+ W^'l^^\xuXi,X2, ...,Xk) 



k+1 
h k-2 



+ E E E wj:iix,,xj)wtr\^i^^K-i) 



m=0 j=l leKj 

r{h) 



d W^fc"\(a;2, . . . , Xj, . . . , Xfe) - wl!'\{x2, . . . ,xi, . . . ,Xk) 



J =2 '^^^ 



V' {xi)Wl{'\xi, ...,Xk)- Ul''\xi; X2, ...,Xk) 



(5.1) 
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m=0 
m=0 j=l leKj 



^ ^t-li^K)-Wl'!^,{xuXK^{j}) (fe) 

(5.2) 



In particular for k = 1. 15.21 reads : 

h-l 

M{x,),M^)Wi''\xi) = Y,wi''-"'\xi)wi'''\x,) + wi''~'\x,,x,) + ul''\x,) (5.3) 

m=l 

It is easy to prove, by double recursion on k and h, that 

k 

Ft\xi, ...,Xk) = 2'' Wt\xi, ,Xk)l[ (5.4) 

1=1 

is a rational function. Introduce the Euclidean division of the polynomial uj.^\xi; xk) 
by M(xi): 

Ul^\xi; xk) = Pt\xi; xk) M{x,) + Qf\x,; xk) (5.5) 

where degP^^'^'* = s — 2 and degQ^.'*'' < d — s. 15.21 becomes: 

V^cr(xi) W^^''\xi,xk) pI:''\xi;xk) 

^/cr{xK) 

_^^'sA wt"'\x,)wl^-\x„XK) wt;;'\x,,x,,XK) 

Mix,) Mix,) 

j!,^ wjf,ix„xj)wtr\^i^^K-i) 

^l^l^l^ M(xi) 

m=0 j=l leKj ^ ^' 

, V- -9 Wf\ixK)-Wf\ix,,XK-{,^ (fe) 

(5.6) 

where the LHS is a rational function of X\ with poles only at the branch points, there- 
fore: 

\fo[x^Wf'' ixx.XK) -^==Pf\xx\XK) 

^aixK) 
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dx 

Res 

x^xi X — Xi 
2s 



.M^)Wt'\x,XK)-Pt\x;xK)) 



= y2 Res^^ (y^wt\x,XK) - Pt\x;xK) 

^ — ' x^ai Xt — X \ 
1=1 ^ 
2s 

= V Res 

i=l 
2s 

+ Res 



dx 



,^j2^t'^\x)wt\x,XK) , w!,'1-'\x,x,xk) 

m=0 



+ 



x^ai Xi — X \ ^ — ' Mix) Mix 

^ m=0 ^ ' ^ 

dx / ^ g ^ wjf,{x ,xj)wtr\^^^K-i) 



i=l 
2s 



Xi — X 



m=0 j=l leKj 



M{x) 



V Res 

' ^ z— >a, 



dx 



i=l 



(5.7) 



xi — X \ ^ — ' dxj 

\ j=2 J 



^ d wl!'\{xK)-Wi''\{x,XK-i,}) ^ Qf\x; xk) 



{xj — x) M{x) 



M{x) 



Two terms in the last line have no pole at Xi = ai, and the other term in the last line 
combines with other lines so as to transform W2 in F2. Thus we get: 



^ 2s dx ^ 
Ft\xi, Xk) = - Res ^^(1 - (5m,o^j,o - ^m.h^i.fc-i) 

2 ^ x^ai Xi — X 



E 



i=\ m=0 J =0 

F^^lix, xi)F^'^~"'\x, xk-i) 



M{x)a{x) 



^ 2s 

+- Res 

2 ^ — ^ x^ai 



dx F, 



k+l 



. xi — X M{x)a{x) 
+Pt''\xi;XK) 



(5. 



where pI:''\xi\xk) is obtained from IH.llH in a way very similar to what we have done 
to leading order. Finally we find: 



Gf\xi,XK) 



rlQ.(^, ^^ "^i+i \x,x,xk) 

s^h-l ( ^ g'^^-"'\x)G^^\x,xk) 



J2i=l x^ai UDi[Xi, X) j^(^), 

+ 2 J2i=l R-6S x^a, J2m=0 dSiiXi, X) — y(^x)dx 

+ Res v-fc-2v- 



dSi{Xi, 



X 



Eh s--\k—2 sr^ 
m=0 2^j=l l^lGKj 

y{x)dx 



(5.9) 

where one should notice that the first line correspond to j = and j = — 1 in the 
second line. 

Let us represent the order A^^^'' A;-loop correlation function G^jl^\xi, . . . ,Xk) as a 
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Figure 1: The 3 skeleton trees contributing to 7^ , i.e. with k + 2h — 2 = A edges. 



black disk with k legs and h holes: 




(5.10) 



Using the Feynmann rules introduced above. can be represented as: 



X 




which means that is obtained by summing over all Feynmann graphs with k 
external legs and h loops. 

The set T^'^^ of all possible graphs with k external legs, and with h loops, can be 
described as follows: First, draw all rooted skeleton trees (i.e. trees whose vertices 
have valence 1,2 or 3), containing fc + 2/i — 2 edges. Draw arrows on the edges, oriented 
from root toward leaves (see figQ). Then draw, in all possible ways, k — 1 external legs, 
and h inner edges, with the constraint that all the vertices of the whole graph have 
valence 3, and so that an inner edge can be drawn only between a vertex and one of 
its descendents (inner edges can never connect different branches of the tree), see fig El 
for the example k = 2, h = 2. Then, each such graph has a symmetry factor. 

We have (see appendix [7|) : 

iVf) := Card rf) = s, (A; -1)! 4^-1 (^^^~^) (^-H) 

where Sh = n['^^ is the number of one-leg graphs in a usual 0^ field theory. The 
generating function s{x) = '^hLi^hx'^'^ is computed in appendix [7| in terms of Airy 



24 




2 ► ' ► I ► 



4 ► ' ► I ► 



■o 



4— H— H— H— H 4 ► ^ ► ^ » 




4 ► ^ ► I ► 



4— H— ► 



Figure 2: All the possible ways of drawing k — 1 = 1 external leg and h = 2 inner 
edges, so that the graphs are trivalent, and that inner edges never connect different 
branches. Notice that all but one graph have symmetry factor 4, and one has 2. 



function. We have: 

Si = 1 , S2 = 5 , S3 = 60 
In particular for genus h = 1: 

Ni^^ :=CardT,^'^ =4'=^^(A;-i; 

Similarly to QUI 



Ten 



(h) 



where the weight W was defined in 14.391 



(5.12) 



(5.13) 



(5.14) 



5.1 Example: One- loop function, genus one 

Let us carry out explicitly the case /c = 1, /i = 1, and recover the result of (THl IT71 IT3] : 
In that case. 15. ll reads: 



M [xi) 



(5.15) 



The RHS is clearly a rational function of xi, and from l2.lil we know that the LHS has 
poles only at the branch-points , and at oo. Introduce the Euclidean division of the 
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polynomial u[^\xi) by M{xi 



U^^\x,)=Pl"{xr)M{x,) + Q'l>{x,) 



(1)/ 



(5.16) 



where degP^^ = s — 2 and degQi^'' < d — s. 



We may thus write: 

^^)w['\xi) - P['\x,) 



Res 



dx 



X — Xi 



{^)Wt'\x)-Pi'\x)) 



dx 



2s 
i=l 

y^Res^^ ^2{x,x)+Q'(\x) 

^ x^ai Xi — 
1=1 
2s 



Xi — X M{x) 
dx W2{x,x) 



(5.17) 



^-^ x^ai Xi — X Mix) 
i=l 



It clearly gives: 



2s 



Gofx, x] 



(5.18) 



Res , , , 

^—^ x^ai y[x)dX 



dSi{xi, x) 



Diagrammatically we have: 





One can check that this result is identical to the function wl'^\x) computed in [TH] . 
Similarly we have: 





1 

^ -2 



+ 2 
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"2 

-or 



=4 ^ ► I ► 




4 



+ 4 -►^ 




+ 4 



and at genus 2 we have: 



~9 



and so on... 



+4 




4 






6 Example, One-cut case s=l, i.e. genus zero curve 

We write: 

a{x) = {x — a){x — h) (6.1) 

It is convenient to map the genus zero hypereUiptical surface into the complex plane 
with the rational map: 

2^(A) = ^ + T(A + A-') (6.2) 

7=^ (6.3) 
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The x-physical sheet is sent to the exterior of the unique disc in the A-plane, and the 
x-second sheet is sent to the interior of the unique disc in the A-plane. We have: 



v/^RA)) = 7(A - A-i) (6.4) 
With this parameterization, all correlation functions are rational functions of the A's. 

6.1 Recursion relations 

All Pfc's are identically vanishing. We have the formula: 



Q ^Jc7{x2) _ 2x1X2- (g+fe) {xi+X2)+2ab 

9X2 (X1~X2) {X1~X2)'^ 




Fk{xi,...,Xk) = -^Res >^>^^- . , . , . \dxfoYk>3 

[x — Xi)M{x)a[x) 



(6.5) 

and ioT k + h > 1. 15.91 gives: 



-Tfc K^I^^K) — l^i=l ^^'^ x^ai {xi-x)M(x)a(x) "'^ 

(m)/ \ j-,{h — m) , 



"^Z^i=l x^ai Z^rn=0 l^j=l l^l€Kj {xi-x)M{x)a(x) ""^ 

+ T^' Res ^fe+i i^^^^^K) ^ 

^ l^i=\ ^^^x^ai (xx-x)M(x)a(x) ""^ 



(6.6) 



6.2 2 point function 

The 2-point function can be written as: 

W^20^i,X2 = -^^In Ai - A2I = 7^7^ In 6.7 

0x1 0x2 0x1 0x2 \Xi — X2 J 

where 

xi=7(Ai + Ar') , X2 = 7(A2 + A2-') (6.8) 



or: 



W2{Xi,X2) 



^ _ ( Va(3:i) - ^Ja{x2) 
^^Jcr{Xl)^Ja{x2) ^ \ Xi - X2 

1 2x1X2 — (a + b){xi + X2) + lab 

2(xi-X2)2 4(xi -X2)V^(2^i)V^(3^2) 



(6.9) 
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In particular we have: 

, (6.10) 
The polynomial A{xi,X2) introduced in 13.311 vanishes identically, and we have: 

W2(x,x) = (6.11) 
All this is well known, see for instance P]. 
6.3 Other correlation functions 

We just give some examples of applications of the general theory described above: 



/ ^ T) f F2{X,X2)F2{X,X3) . 

Fslxi, X2, xs) = - Res ^ , . dx 

a,b \[x — xi)M[x)a[x) 

F2{a,X2)F2{a,X3.) 



(6.12) 



(a — xi)M{a){a — h) 

F2{b,X2)F2{b,X3) 

{b-xi)M{b){b-a) 
b — a 

(a — xi){a — X2){a — X3)M{a) 
b — a 

~ {b - xi){b ~ X2){b - x^)M{b) 



I.e. 



vv^[xi, X2, X3) - [0- a) „ / , X / / X r-r =^ {b.i6) 

8^^a{Xl)^^(T{X2)^/c^{X3) 

which is the usual of pp. 

(6 — a)^ dx 



Res 



32 a,b (xi - x)M{x)a{x) 



2 



(b — aY ^ dx 1 
Res 



32 a (x - a)2 {xi - x)M(x) (x - b) 

(b — a)^ ^ dx 1 
+ ^ ^ Res 



2 



32 6 (x-6)2 (xi -x)M(x)(x-a) 
(6-a)2 / 1 



2 



32 \{xi - x)M{x){x -by 
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(6-a)2 / 1 



) 



32 \{xi- x)M{x){x - ay 
{h-af ( 1 -2xi -6 + 3a 



M\a) 1 



x=b 



) 



32 \M{a) (a-6)3(xi -a)2 
(6 - a)2 / 1 -2x1 - a + 36 



M(a)2 (a-6)2(xi - a) 
M'ib) 1 



) 



32 V^(6) {h-af{xi-hY 



MipY {h-af{xi-h) 



(6.14) 



which again agrees with and other results in the hterature. 

7 Conclusions and prospects 

In this article, we have found a 0^ Feynmann graph formulation for computing all 
correlations functions to all powers of in the one-hermitian matrix model. First, it 
would be interesting to find out to which field theory it corresponds. One is tempted 
to compare with Liouville's theory (which is not cubic) or to a fermionic theory. 

We claim that this approach is more efficient for actual calculations, than the 
method existing previously in the literature jSISl- Indeed, in [20], one has to construct 
the correlation functions recursively, by expanding them on basis functions which are 
themselves constructed recursively by taking derivatives with respect to the potential. 
For instance, one does not get any simplification in the method of |21 El by assuming 
an even potential, or by assuming a quadratic potential. The method presented here, 
works for fixed potential (for instance quadratic), and does not need to construct any 
basis of functions. 

Another important point for the method presented here, is that it is expressed 
in terms of geometrical fundamental objects on the spectral curve. This is another 
evidence of the deep link between tau functions and complex geometry. 

There are other expressions in the literature involving Residues of geometrical ob- 
jects (for instance [331211111 0113 113 I^Hllini), namely, only the Bergmann kernel and 
not the abelian differential. However, we claim that it should be simpler to compute 
the residue of a function with a simple pole (the abelian differential), than the residue 
of a function with a double pole (Bergmann kernel). 

Moreover, the whole procedure described here, can be applied with very small 
adaptations to other matrix models, in particular the 2-matrix model, and to non- 
hermitian matrix models (in particular /5 = 1, 2, 4 models), this work is in progress and 
will be available shortly [201 • In the 2-matrix model with potentials of degree di + 1 
and (i2 + 1, the computation of correlation functions of the first matrix only involves 
(^2 vertices (i.e. cubic, quartic, d2 -|- 2-legs-vertex) , instead of only one cubic vertex 
equal to l/2y{x) for the 1-matrix model. This will be further explained in |2Uj . 
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The observable we have not computed in this article is the free energy: 

^ oo 

G,:=-—\r.Z:=Y,N-'''Gt^ (7.1) 
The free energy does not appear in the loop equations. It satisfies: 

oV[xi) 

therefore, in order to compute the free energy, one has to integrate with respect to the 
potential, i.e. one can no longer keep the potential constant. One would reasonably 
make the following conjecture for h > 2: 

h;(t) (7.3) 



for example for h = 2 






+ 



unfortunately, these Feynmann graph don't make sense (the abelian differential diverges 
at coinciding points). The conjecture is that the G'^q^ are related to traces of powers 
of the laplacian on the spectral curve. For instance it is known that G^q^ is related to 
the determinant of the laplacian |15j . 
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Appendix A Cardinal of Tj^^^ 

The cardinal of 7^ and of T^'^^ can be computed by setting W = 1 in 14.401 and in l5.14[ 
and then using the recursion relations 14.351 and 15.91 

One thus gets for /c > 3: 

No = l 



7=1 ^ ^ 



31 



writing: 

ro:=0 , ri := 1 , := — ^ 

lA. II becomes for k > 2: 

k 

i=o 

We introduce the generating function: 

oo 

R{x) ■=^rkx'' 

k=0 

and thus IA.3I becomes: 

R{x) - x = R^{x) 

whose solution is: 



k=i 



which imphes: 



r, = -t^ 1)^+1 2^^-^ ^ ^"^^ • • • (I - ^) _ - 3)! 



2 V^/ 
and thus, we obtain lOTI 

Arfc = 2^-2(2A;-5)!! 



A; -2! 



For higher genus, we have for k > 1 and /i > 1: 

j=0 m=0 V / 

writing: 
]9l becomes for k > 0, h> 1: 



j=0 m=0 



We introduce the generating function: 

oo 

Rh{x) := 5^ 



W k 
rl X 



k=0 

32 



and thus lA. lH becomes for h > 0: 

h 

Mx) = R{x) , Rh{x) = R'f^_,{x) + Rm{x) Rh-m{x) (A.13) 

m=0 

which can also be written for /i > 1: 

h-l 

(1 - 2R{x)) Rh{x) = + ^ R^{x) RH-m{x) (A.14) 

m=l 

using \KM. it is easy to see, by induction on h that for h>l one has: 

Rh{x) = Sh (l-4x)"— (A.15) 
where the coefficients Sh obey for /i > 1: 

Si = l , Sh = 2{3h- A)Sh-l + YsmSh-m (A.16) 

m=l 

or, if we define Sq := — |, it can be written for any h > 1: 

h 

= 2{3h-4)Sh-l + J2^rnSh-m (A.17) 

m=0 

we introduce the generating function: 

oo 

S{x) ■.= J2^hx'' (A.18) 



h=0 



it obeys: 
If one writes 
and 

one has: 

whose solution is; 



= S\x)-^ + 6x^S'{x)-2xS{x) (A.19) 



X' 



-2/3 



e = ^ (A.20) 



S{x) = -xh{^) (A.21) 
^ = h\0 + h\0 (A.22) 



(A.23) 
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or 



and thus: 

S(x) = — ^ . = — 1 + ^ ^ A.24 

^/ V 1 ^ f tdtc 2 c s / » 

^(^) = "2 - ^ ^ (A-25) 
J ate 2 e s 

S{x) = --+ ^'=" -^(^ +^)' (A.26) 

yr^oo (6fc+3)!! v^oo j;'' (li/i+l)! 

c,/ I 2^/1=0 32^+1 (2h+l)! _ 1 Z^/i=0 S^h+i 23''+2 (2/1+1)! (3/1+2)! A oyN 

W ~ ~2 1 _L xfe(6fe-l)!! ~ ~^ ^^^^ (6/t)! i^-^'J 

+ Z^h=l 32ft (2/i)! Z^/i=0 32h 23ft (2/i)!(3/i)! 

In the end we have: 

/ 3(/t~ 

flarH T:^'''' = Ci- - 1 i'^-l I 2 

A;- 1 

The whole function is thus: 

k h h 



TVf ^ := Cardrf ) = s/, (A; - 1)! 4^-^ f '^'2 " + f M (A.28) 



^2 



e 3-64 



t3 _i(i-4x) 



(A.29) 
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